arXiv: 1501.07367v3 [cond-mat.str-el] llJun2015 


Journal of the Physical Society of Japan 


FULL PAPERS 


Electron Correlation Effects in Non-Centrosymmetric Metals 
in the Weak Coupling Regime 

Daisuke Maruyama* and Youichi Yanase'’^’^ * 

^Graduate School of Science and Technology, Niigata University, Niigata 950-2181, Japan 
^Department of Physics, Niigata University, Niigata 950-2181, Japan 
^Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan 

The two-dimensional Rashba-Hubbard model is investigated in order to clarify the electron correlation effects in non- 
centrosymmetric metals. The renormalization effect on Rashba spin-orbit coupling (RSOC) is calculated on the basis of 
second-order and third-order perturbation theories. We show that RSOC is enhanced by the electron correlation. On the 
other hand, the spin-splitting of the Fermi momentum (SFM) is robust against the electron correlation effect because the 
enhancement of RSOC is cancelled by the k-mass renormalization. The cancellation is almost perfect in often adopted 
models in which the momentum dependence of RSOC is linear in the quasiparticle velocity. A small correction to the 
SFM appears otherwise. We show that the same results are obtained for other antisymmetric spin-orbit couplings in 
non-centrosymmetric systems, such as the synthetic spin-orbit coupling in cold atoms. 


1. Introduction 

A breakthrough in exploring exotic states of matter in¬ 
duced by spin-orbit couplings has been one of the high¬ 
lights of recent condensed matter physics. In particular, non- 
centrosymmetric metals lacking an inversion symmetry in 
the crystal structure exhibit various intriguing phenomena 
such as spintronics,^’^ chiral and helical magnetism,^’and 
non-centrosymmetric superconductivity.^ An antisymmetric 
spin-orbit coupling relating the momentum and spin of elec¬ 
trons appears owing to the lack of inversion symmetry, and 
it plays an essential role in these phenomena.^ Although 
various antisymmetric spin-orbit couplings exist depending 
on the symmetry of crystals,® Rashba spin-orbit coupling 
(RSOC)^ has been investigated most intensively, probably be¬ 
cause it appears not only in the bulk® " but also in artificial 
heterostructures.The antisymmetric spin-orbit coupling 
gives rise to the spin-splitting of Fermi surfaces, and thus, a 
single-particle state acquires a spin texture in the momentum 
space.® Such an electronic structure results in exotic quan¬ 
tum phases and quantum transport. Spin-split Fermi surfaces 
have been observed by angle-resolved photoemission spec¬ 
troscopy (ARPES)"’and by de Flaas-van Alphen (dHvA) 
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measurement. 

The interplay between the antisymmetric spin-orbit cou¬ 
pling and electron correlation effects often plays an im¬ 
portant role. For instance, unconventional superconductivity 
with a giant upper critical field occurs in strongly corre¬ 
lated non-centrosymmetric systems.Furthermore, it has 
been shown that novel electromagnetic responses, such as 
a magnetoelectric effect and anomalous Hall effect, are en¬ 
hanced by the electron correlation.^^ Thus, the interplay be¬ 
tween the antisymmetric spin-orbit coupling and electron cor¬ 
relation appears to be an important issue. From the theo¬ 
retical point of view, unconventional long-range order, such 
as magnetism^®^^® and superconductivity,^^’^®’^®^®® has been 
intensively clarified. Furthermore, spontaneous inversion- 
symmetry breaking and the emergence of an antisymmetric 
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spin-orbit coupling due to electron correlation effects have 
been studied.®'®® However, quasiparticles renormalized by 
electron correlation effects in the Fermi liquid state have only 
been investigated in a few works.®*’®®’®" 

A Fermi liquid theory for non-centrosymmetric metals 
(chiral Fermi liquid theory) has been formulated on the ba¬ 
sis of a two-dimensional Rashba-Hubbard model by Fuji- 
moto,®® although only the diagonal self-energy was explic¬ 
itly calculated on the basis of second-order perturbation the¬ 
ory. Although a phenomenological spin fluctuation model in 
three dimensions has been investigated by Tada et al.,®*’®" 
again the off-diagonal self-energy was neglected. As shown 
in Ref. 22, the renormalization of spin-orbit coupling is ne¬ 
glected in these calculations. A one-dimensional Rashba- 
Hubbard model has been analyzed by Goth and Assaad,®® 
but Landau quasiparticles break down in one-dimensional 
systems. In this paper, we clarify electron correlation ef¬ 
fects in the chiral Fermi liquid state on the basis of the two- 
dimensional Rashba-Hubbard model by calculating both di¬ 
agonal and off-diagonal self-energies on an equal footing. The 
calculation relies on the perturbation expansion with respect 
to the Coulomb interaction, and therefore, we obtain reliable 
results in the weak coupling regime. 

Two-dimensional electron gases in semiconductors have 
been theoretically studied for a similar purpose.®® Also, quasi¬ 
particle properties in the presence of antisymmetric spin- 
orbit coupling and a screened long-range Coulomb interaction 
have been investigated. It has been shown that the spin-orbit 
coupling is enhanced by the momentum dependence of the 
screened Coulomb interaction,®" but the quasiparticle proper¬ 
ties are hardly affected by the interplay between them.^"^® 
Interestingly we obtain similar results for correlated met¬ 
als, although many-body effects arising from the short-range 
Coulomb interaction are the main subject of the paper. We 
also show the important role of anisotropy in the electron dis¬ 
persion relation, which is not negligible in metals. 

The paper is organized as follows. In Sect. 2, we intro¬ 
duce the Rashba-Hubbard model and formulate the renormal¬ 
ization of quasiparticles in the presence of electron correla- 
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tion and spin-orbit coupling. Numerical and analytic results 
of second-order perturbation theory are shown in Sect. 3. In 
Sect. 4, we show the results of third-order perturbation theory 
in order to clarify higher-order corrections. A brief summary 
and discussion are given in Sect. 5. 

2. Formulation 

In this section, we introduce the Rashba-Hubbard model 
and formulate the Green function, self-energy, and effective 
mass of quasiparticles in the presence of RSOC. 

2.1 Rashba-Hubbard model 

The effects of electron correlation in metals lacking an 
inversion symmetry are studied on the basis of the two- 
dimensional Rashba-Hubbard model, 

( 1 ) 

■(Tc\cks', (2) 

k,s k,s.s' 

f^int = (3) 

i 

where Cks (clj is the annihilation (creation) operator for an 
electron with momentum k and spin s i- The electron 
number operator for spin s at site i is denoted as 

We consider a simple square lattice and assume a tight- 
binding model, e(k) - -2fi(cos k^-Hcos ky)-i-4t2 cos cos ky- 
fj., by taking account of the nearest-neighbour and next- 
nearest-neighbour hoppings. The chemical potential ju is in¬ 
volved in the dispersion relation. The second term in the 
single-particle part Hq represents the antisymmetric spin- 
orbit coupling arising from the lack of inversion symme¬ 
try. Following the conventional notation,^ it is characterized 
by the g-vector g(k). We focus on RSOC, which has been 
studied in various fields of condensed matter physics. Thus, 
in Sects. 3 and 4 we will assume a g-vector that repre¬ 
sents RSOC. Although the momentum dependence of the 
g-vector is determined by the orbital wave function in the 
Bloch state,^^ we will adopt a simple form of RSOC, g(k) - 
2fi(- sinky, sinker, 0), which is justified in the absence of or¬ 
bital degeneracy in the electronic structure.On the other 
hand, we will discuss other kinds of antisymmetric spin-orbit 
coupling in Sect. 5. Thus, we adopt a general form of the g- 
vector in this section. 

The interacting part Hint represents the on-site Coulomb 
repulsion. We investigate electron correlation effects in the 
weak coupling regime by using the perturbation expansion 
with respect to the Coulomb interaction U. 


where I is the 2x2 unit matrix and a>„ = {2n - 1)7 tT is 
the fermion Matsubara frequency. The noninteracting Green 
function in the chirality basis is obtained as 


Gf{k,icJn) 


1 

\Un - E{k) - Aa\g{k)\' 


(5) 


with A - ± being the chirality index. 

The dressed Green function is obtained by taking account 
of the self-energy £(k) = [2js-(k)]. As we will show later, not 
only the diagonal self-energy Y.a-a-{k) but also the off-diagonal 
self-energy ’Eg-a-ik) plays an important role, although the latter 
was neglected in previous works.Carrying out analytic 
continuation, the retarded Green function is described as 


G\k) 


I M-e-ag,-X’L -a(gx - i.gv) - \ 

^ -a(g, + igy) co-e + ag,- j ’ 


where k - (k, ai). Adopting the vector representation of the 
self-energy. 


rR 


(7) 


where + E^j^)/2, E^ = + 
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E^p/2, E^ = (E« - Ej5p/2i, and Ef = (E^^^ - E«p/2, the 
dressed Green function is represented in a familiar form. 


G^(k) 


Ico-e'-ag'^ -a(g'^ - ig'y) \ 
\ -o'ig'x + i^y) oj-e' -hag'J 


Z 

A=± 




(T 


G^(k). 


( 8 ) 


The renormalized Green function in the chirality basis is ob¬ 
tained as 




1 

u - s’{k) - Aa\g'{k )\’ 


(9) 


where e'{k) = e(k) H- E^ (k). The renormalization of spin-orbit 
coupling is taken into account by the self-energy correction of 
the g-vector, ag'{k) = ag{k) -H ReE*^(k)."^^ Here we dropped 
the imaginary part of the spin-dependent self-energy ImE’^(k), 
because it is negligible at low temperatures in a Fermi liquid 
state, ImE’^(k) oc T^. Figure 1 shows the self-energy obtained 
by second-order perturbation theory, and we indeed see the 
negligible imaginary part around ai - 0. This property is not 
altered in third-order perturbation theory, indicating the Fermi 
liquid state. 

The Fermi surfaces of spin-split bands are defined by the 
singularity of the Green function and are thus obtained by 
solving the equation 


e'(A:F^,0) + da|g'(A:F^, 0)1-0. (10) 


2.2 Green function 

In this subsection, we formulate the renormalization of 
quasiparticles due to the electron correlation effect. In the 
presence of spin-orbit coupling, the noninteracting Green 
function is described in the matrix representation as 


G<°>(jt,iw„) 


hoJn-B- ag, -aigx - igy) \ ' 
\ -a(8x + igy) ioJn-s + ag^ j 


-Z 


7 + ^ri 

lgl 


cr 


Gf{k,iajn), 


(4) 


The Fermi momentum of the d-band is denoted as kp^ and 
the spin-splitting of the Fermi momentum (SFM) is defined 
as A^f - I^ff - ^F-l. 

We also take into account the correlation correction of the 
chemical potential. First, we calculate the chemical potential 
fiQ 3.1 U - 0 for which the electron density per site is n. Next, 
we calculate the self-energy, and the chemical potential is cor¬ 
rected as ju = juo H- dfi, for which the dressed Green func¬ 
tion leads to the electron density, lim;,_F+oTr Zk G(k)e'""'' = 
TTZkG(k)+i =n- 


2 











J. Phys. Soc. Jpn. 


FULL PAPERS 



Fig. 1. (Color online) Frequency dependence of retai'ded self-energies at 
temperature T = 0.01 obtained by second-order perturbation theory. We show 
the results at Fermi momenta k = ftF+ II [HO] for {t\J 2 ) - (CO), f/ = 4, 
and O' = 0.1. Solid (dashed) lines show the real (imaginaiy) part of (a) Sq(/:), 
(b) Sj(A:), and (c) (^). We adopted the Fade approximation for the analytic 

continuation. 


tial role in correlated non-centrosymmetric metals as we will 
show below. 

We adopt an approximation formula for a numerical calcu¬ 
lation. We obtain the self-energy at cu = 0 through the Mat- 
subara self-energy 




'La(k, Un - 7 tT) + I.a(k, Wn = -TtT) 

2 


■ (14) 


The temperature is assumed to be T = 0.01 in the follow¬ 
ing numerical results. Frequency derivatives are calculated by 
using the Kramers-Kronig relation as 


5ReS|J(jt,w) 


do) 


dlrc&aik, Wn) 


fj=0 


dUn 


a>n=0 


ImS„(^, Wn = ttT) - ImSQ.(^, Wn = -nT) 


(15) 


We also calculated the retarded self-energy 2^(A:) using the 
Fade approximation (see Fig. 1) and estimated the w-mass 
and A:-mass. It has been conhrmed that the two numerical es¬ 
timations coincide with each other. Thus, the approximation 
formulas Eqs. (14) and (15) have been justihed. 


3. Second-Order Perturbation Theory 

We here show the results of second-order perturbation the¬ 
ory for the two-dimensional Rashba-Hubbard model. Since 
the hrst-order self-energy is involved in the correction of the 
chemical potential, second-order perturbation theory is the 
lowest-order theory justihed in the weak coupling regime. 
We calculate the second-order self-energy represented by the 
skeleton diagrams in Fig. 2. In the absence of RSOC, the self¬ 
energy corrections represented by Fig. 2(2B) disappear. On 
the other hand, these terms give rise to the renormalization 
of RSOC. Thus, we take into account both Figs. 2(2A) and 
2(2B), although the latter was neglected in a previous study. 


2.3 k-mass and io-mass 

Next, we introduce the effective mass of quasiparticles. The 
effective mass of the d-band m*^ is obtained as the product of 
the w-mass and k-mass. 


mx mx mx ’ 


(11) 


where niA is the bare mass of the d-band. The w-mass is given 
by the frequency derivative of the self-energy. 


nq 


dReiq(kp,,co) 


dco 


-Aa 




a >=0 


d(jj 


( 12 ) 


CO-Q 


Furthermore, quasiparticles acquire a k-mass renormalization 
through the momentum derivative of the self-energy. 


m\ 

niA 


__OK_ 

de(k) , (JReElfCt.O) , , „,fl|g'(t,0)| 

— + — M — + 4a 


(13) 


The k-mass renormalization is often neglected in Fermi liq¬ 
uid theory for strongly correlated electron systems because 
it is quantitatively less important than the w-mass renormal¬ 
ization. However, the k-mass renormalization plays an essen- 



Fig. 2. Skeleton diagrams of self-energy in second-order perturbation the¬ 
ory. The solid and wavy lines show the bare Green function G®(k) and 
Coulomb interaction U, respectively. 


3.1 Numerical results 

First, we show the numerical results of the effective mass, 
the renormalization of RSOC, and the SFM. We choose the 
coupling constant of RSOC as a = 0.1 throughout the pa¬ 
per. We assume U - A unless otherwise specihed. In this 
subsection, the next-nearest-neighbour hopping is neglected 
for simplicity, and thus we assume (fi,f 2 ) = (1>0). Then, 
the Rashba-type g-vector is represented by the velocity as 
g{k) - 2fi(-sinkj,,sinki,0) = [-v'y(k:), V;c(k:), 0]. In the next 
subsection (Sect. 3.2), we will show that this relation between 
the RSOC and the velocity plays an essential role. 


3 


























J. Phys. Soc. Jpn. 


FULL PAPERS 


Figure 3 shows the w-mass as a function of the electron 
filling 11 . As is known from Fermi liquid theory, the effective 
mass is enhanced by the electron correlation through the oj- 
mass. The mass enhancement is pronounced near the half¬ 
filling, « = 1, because of the large density of states (DOS) due 
to the Van-Hove singularity at k - (tt, 0) and (0, tt). Because 
we assume a small RSOC compared with the Fermi energy, 
the effect of RSOC on the w-mass is negligible. Indeed, the 
three lines in Fig. 3 almost coincide with each other. Thus, 
the band dependence of the effective mass is small in the weak 
coupling regime. 


2 
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through the spin-dependent part of the self-energy E(k). Be¬ 
cause the main contribution to S(k) comes from the diagram 
represented in Fig. 2(2B), it is essential to take into account 
these terms when studying of correlation effects on the RSOC. 



Fig. 5. (Color online) Filling dependence of the renormalization of RSOC 
obtained by second-order perturbation theory. Squares and circles show 
\g'{k,cij = 0)|/|^(A:)| at the Fenni momentum k = || [110] for /i = + 

and A = respectively. 


Fig. 3. (Color online) Filling dependence of the tu-mass at A: = II [110] 
obtained by second-order perturbation theory. Squares and circles show the 
tu-mass of the /I = + and - band, respectively. The tu-mass in the absence of 
RSOC (n = 0) is shown by crosses. 

In contrast to the oi-mass, the k-mass is suppressed by elec¬ 
tron correlation effects. Indeed, Fig. 4 shows that < 1. 

However, the total effective mass [Eq. (11)] is enhanced be¬ 
cause the w-mass renormalization is much larger than the k- 
mass renormalization. For lal 1, both w-mass and k-mass 
renormalization mainly originate from the diagonal compo¬ 
nent of self-energy E^^^k) represented by Fig. 2(2A). 



Fig. 4. (Color online) Filling dependence of the X:-mass at A: = |! [110] 

obtained by second-order perturbation theory. Squares and circles show the 
A:-mass of the d = + and — band, respectively, while crosses show the A:-mass 
in the absence of RSOC (a = 0). 

Now we discuss the renormalization of RSOC due to the 
electron correlation effect, which is given by \g'{k,u - 
0)|/|g(k)|. Figure 5 shows that \g'{k,cij - 0)|/|g(k:)| > 1 
around the Fermi surface irrespective of the electron den¬ 
sity, and thus RSOC is enhanced by the electron correlation 


For our choice of hopping integrals (fi,f 2 ) = (CO), the 
Hamiltonian has a particle-hole symmetry because e{k) +i-i- 
-e(Q -k)- ji and g{Q -k) — g(k) with Q - (n, n). Then, the 
particle-hole transformation, c*, —> Cg changes the signs 
of the chemical potential (jj —> -/r) and RSOC {a —» -a), 
and thus the electron density is changed as « ^ 2 - n. 
Therefore, the renormalizations of the w-mass, k-mass, and 

RSOC show symmetric behaviors ^ l«=2-«o 

\g'{kY,l,Q)\l\g{kvA)\n=na = \g'{.ky^A,Q)\l\gi.kv-A)\n=2-n„ (sCe 
Figs. 3-5). 

The renormalization of RSOC is anisotropic in the momen¬ 
tum space as shown in Fig. 6. Generally speaking, a large en¬ 
hancement of RSOC occurs in the vicinity of the Fermi sur¬ 
face, particularly at the momentum k = kp || [100] and [010]. 
The spin texture indicated by the direction of the g-vector is 
also changed by the electron correlation effect (not shown). 
Although the correlation effect on the spin texture is negli¬ 
gible in the weak coupling regime, the spin texture may be 
considerably changed in the strong coupling regime. 

Finally, we show the SFM Akp. Although the enhance¬ 
ment of renormalized RSOC implies an increase in Akp, the 
SFM remains unchanged by the electron correlation effect, as 
shown in Fig. 7. Note that the scale of the horizontal axis is 
chosen to be the same as in Fig. 5. We can not observe the 
electron correlation effect on Akp at this scale. Thus, it is indi¬ 
cated that the enhancement of RSOC is compensated by other 
effects. The k-mass renormalization indeed cancels out the en¬ 
hancement of RSOC, as we show in the next subsection. 

3.2 Analytic calculation 

We here clarify the cancellation between the k-mass renor¬ 
malization and the renormalization of RSOC. In this subsec¬ 
tion, we adopt the lowest-order theory with respect to RSOC, 
which is justified for a weak RSOC, a|g^(k)| «: ep, as realized 
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Fig. 6. (Color online) Momentum dependence of the renormalization in 
RSOC. We show \g'{k, cj = 0)|/|g(A:)| for fillings (a) n = 0.4, (b) 0.6, (c) 0.8, 
and (d) 1.0. Solid lines show the spin-split Fermi surfaces. 



Fig. 7. (Color online) {/-dependence of SFM A^f for k || [110]. We 
plot AA:f for various electron densities n =0.2, 0.4, 0.6, 0.8, and 1.0, but 
AA:f(U) ^ AA:F(t/ = 0) irrespective of the electron density. The scale of the 
horizontal axis is chosen to be the same as that in Fig. 5 in order to show the 
almost complete cancellation between the A:-mass renormalization and the 
enhancement of RSOC. 


in most non-centrosymmetric metals. The SFM is obtained as 




2 a|g'(^F, 0 )l 
\vHkFM ’ 


(16) 


where kp is the Fermi momentum at o- = 0 and w) = 

\k=kF is the Fermi velocity renormalized by 
the A:-mass. Thus, the SFM is affected by the renormalization 
of the A:-mass and RSOC but not affected by the ru-mass. 

The perturbation expansion in terms of RSOC is carried out 
by expanding the noninteracting Green function as 

Gf G^”^ + G<°^da|g|G<°’ + G®Ta|g|G('’^Ta|g|G('’^ + • • • , 


(17) 


where G^^\k) - {u) - £(^)]“'. Up to the first order in a, the 
diagonal self-energy So(A:) is obtained from the diagram in 
Fig. 2(2A). It is the zeroth-order term with respect to a. On 
the other hand. Fig. 2(2B) represents the first-order terms in 
the off-diagonal self-energy Ej:(A:) and 2v(A:). Thus, the lowest- 


order terms of the self-energy are obtained within the first or¬ 
der of a. 

Differentiating the diagonal self-energy, we obtain the 
renormalized Fermi velocity as 

y|^(^) = vAk) + t/2Re J] 0(^)G®(^ - kfvAq - k), (18) 

v^(^) = Vy{k) + G^Re ^ - k)\(q - k), (19) 

where 

0 (^) = ^G®(U)G<°)(^-U), ( 20 ) 

k’ 

and v„(^) = de{k)ldka. On the other hand, the renormalized 
g-vector is obtained as 

g'Ak) = 8x(k) + G^Re ^ cf>(q)G^°\q - kfgAq - k), (21) 

g'Ak) = gy{k) + G^Re y] 0(^)G®(^ - kfgy(q - k). (22) 
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The derivation of Eqs. (18)-(22) is given in the Appendix. 
Now the similarity between the Fermi velocity and RSOC is 
clear. The Fermi velocity and RSOC acquire the same renor¬ 
malization when the relation g{k) - C{—Vy{k), Vx{k), 0] is 
satisfied with C being an arbitrary constant. This relation 
is indeed satisfied in the Rashba-Hubbard model adopted in 
Sect. 3.1. Then, the SFM is not renormalized by the electron 
correlation effect. 


2a|g'(A:F,0)| ^ 2a|g(A:F)| 
|vH*^f,0)| |v(A:f)| 


(23) 


The cancellation is not complete owing to higher-order terms 
with respect to the RSOC a, but we see almost complete 
cancellation between the A:-mass renormalization and the en¬ 
hanced RSOC for a moderate RSOC of a = 0.1 (see Fig. 7). 



Fig. 8. (Color online) SFM Akp for n =0.6, 0.4, 0.2, 0.8, 1.0, 1.8, 1.6, 
1.2, and 1.4 from top to bottom. We assume (fi,/ 2 ) = (1,0.3), although we 
assumed (/], ^ 2 ) = (1,0) in Fig. 7. The other parameters are the same as those 
in Fig. 7. 


The above results imply that the SFM is renormal¬ 
ized by the electron correlation when the relation g{k) - 
C[-Vy(k), Vx(k), 0] is not satisfied. For instance, we can 
choose the parameters (fi,f 2 ) = (C0.3) so that g{k) 96 
[-Vv(^), Vj:(^), 0]. Indeed, Fig. 8 shows that the SFM is af¬ 
fected by the Coulomb interaction. A finite correction to the 
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SFM is obtained, although it is substantially reduced by the 
cancellation of the A:-mass renormalization and the enhance¬ 
ment of RSOC. The sign of the correction depends on the 
band structure. We see that the SFM is enhanced (suppressed) 
by the electron correlation effect when the Fermi surface is 
electron-like (hole-like). This particle-hole asymmetry in the 
correction to the SFM is caused by the next-nearest-neighbour 
hopping t 2 , which induces the particle-hole asymmetry in the 
band structure. However, we note that the robust SFM inde¬ 
pendent of U (see Fig. 7 for example) does not require the 
particle-hole symmetry in the band structure. As we showed 
above, the SFM is not renormalized by the electron correla¬ 
tion when g{k) - C[-Vv(k), Vi(k), 0]. We can choose the g- 
vector so as to satisfy this relation even when t 2 + 0 . 

4. Third-Order Perturbation Theory 

So far we have investigated the weak coupling region of the 
Rashba-Hubbard model on the basis of second-order pertur¬ 
bation theory. In this section, we examine higher-order cor¬ 
rections by comparing third-order perturbation theory with 
second-order perturbation theory. The third-order terms of 
self-energy are diagrammatically represented in Fig. 9. These 
terms are classified according to the leading order with re¬ 
spect to the RSOC a. The zeroth-order terms of a are 
and while and £(3F)-(3J) are first-order terms. Since 
E(3D) £( 3 E) ^j.g higher-order terms, they are negligible for 

Q- = 0.1. 



Fig. 10. (Color online) {/-dependence of the ^J-mass in the A = + band 
at k = kp-^ II [110]. Crosses, stars, and circles show the results of third- 
order perturbation theory for n =1.0, 0.8, and 0.6, respectively. The results 
of second-order perturbation theory are shown for n = 1.0 (squares), 0.8 
(triangles), and 0.6 (diamonds) for comparison. 





(3H) 






Fig. 9. Diagrammatic representation of third-order terms of self-energy. 


First, we assume the dispersion relation (fi, f 2 ) = (1,0), as 
in Sect. 3.1, and calculate the renormalization of the effective 
mass, RSOC, and SFM. Figures 10 and 11 show the w-mass 
and k-mass on the A = + band, respectively. As shown for 
the Hubbard model without spin-orbit coupling,^® the third- 
order correction partly cancels the second-order terms. In¬ 
deed, both w-mass and k-mass renormalization are suppressed 
by the third-order terms. The third-order terms give rise to a 
particularly large correction to the k-mass. A special case ap¬ 
pears at half-filling, ii - 1. In this case, the third-order terms 
are negligible because the two leading-order terms and 
E*^®' cancel each other because of the particle-hole symme¬ 
try."*® 

The renormalization of RSOC on the Fermi surface of the 
A = + band is shown in Fig. 12. It is shown that the third- 


Fig. 11. (Color online) {/-dependence of the /:-mass in the /I = + band at 
k = kp+ II [110]. Symbols indicate the same electron density and the same 
order of perturbation theory as those in Fig. 10. 


order correction reduces the enhancement of RSOC. Inter¬ 
estingly, the SFM is invariant against the electron correla¬ 
tion even when we take into account third-order terms (see 
Fig. 13). Thus, the cancellation between the renormalization 
of RSOC and the k-mass is not an artifact of second-order 
perturbation theory. Indeed, the cancellation occurs in each 
order of U when g(k) - C[-Vv(k:), Vji(k:), 0]. Therefore, it is 
expected that the SFM is robust against electron correlations 
not only in the weak coupling regime but also in the strong 
coupling regime. 

Next, we choose the parameters (fi,f 2 ) = (1,0.3), as in 
Fig. 8 , so that g(k) 9 ^ [-Vj,(k:), V;c(k:), 0]. Then, the SFM is 
renormalized by the electron correlation effect as shown in 
Fig. 14. Thus, we obtain qualitatively the same results as 
those in second-order perturbation theory, although the elec¬ 
tron correlation effects are reduced by the third-order correc¬ 
tion terms. 

5. Summary and Discussion 

Exotic quantum phases and intriguing electromagnetic re¬ 
sponses are induced by the antisymmetric spin-orbit cou¬ 
pling in non-centrosymmetric metals.'® Naturally, the in¬ 
terplay between electron correlation effects and spin-orbit 
coupling in non-centrosymmetric metals has attracted inter¬ 
est. We have investigated electron correlation effects in the 


6 



























J. Phys. Soc. Jpn. 


FULL PAPERS 



Fig. 12. (Color online) t/-dependence of the renoiinalization of RSOC at 
k = A:f+ II [110]. Symbols indicate the same electron density and the same 
order of perturbation theory as those in Fig. 10. 
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Fig. 13. (Color online) SFM obtained by third-order perturbation theory 
for(ri,r 2 ) = (1,0). We show the results for electron densities « = 0.2,0.4,0.6, 
0.8, and 1.0. The momentum is chosen to be parallel to the [110] axis. The 
SFM is invariant against the electron correlation effect on the whole Fermi 
surface irrespective of the electron density. The scale of the horizontal axis is 
chosen to be the same as that in Fig. 12. 



Fig. 14. (Color online) SFM obtained by third-order perturbation theory 
for = (1,0.3). for« =0.6, 0.4, 0.8, 0.2, 1.8, 1.0, 1.6, 1.2, and 1.4 
are shown from top to bottom at f/ = 2. 


two-dimensional Rashba-Hubbard model, wbicb is a minimal 
model for a non-centrosymmetric (cbiral) Fermi liquid. Tbe 
effective mass, tbe renormalization of spin-orbit coupling, and 
tbe spin-split Fermi surfaces were calculated on tbe basis of 
perturbation theory. We showed that the electron correlation 
enhances the spin-orbit coupling but the SFM Ak-p is almost 


invariant against Coulomb interaction. 

In second-order perturbation theory, both numerical and 
analytic calculations show that the enhancement of spin-orbit 
coupling is cancelled by the A:-mass renormalization. Thus, 
Akp is not renormalized by the electron correlation. The can¬ 
cellation is complete when the g-vector of RSOC is repre¬ 
sented by the velocity as g{k) - C{—v^{k),vJJ<.),Qi\. Other¬ 
wise, a finite correction to the SFM appears, but it is reduced 
by an incomplete cancellation. 

We numerically examined the validity of second-order per¬ 
turbation theory by calculating third-order correction terms 
with respect to the Coulomb interaction. Generally speak¬ 
ing, third-order terms partly cancel the leading-order second- 
order terms.Indeed, the renormalization of the effective 
mass and spin-orbit coupling is decreased by the third-order 
terms. The SFM remains invariant against the electron cor¬ 
relation in third-order perturbation theory when the relation 
g{k) - C{—Vy{k), vjji), 0] is satisfied. Thus, the robustness 
of the SFM against the electron correlation is not an artifact 
of second-order perturbation theory and is expected to be an 
exact property. 

Although we considered a Rashba-type spin-orbit coupling, 
our results are generally valid for other kinds of antisym¬ 
metric spin-orbit coupling. For example, we confirmed that 
the spin-orbit coupling is enhanced but the SFM is invari¬ 
ant against the electron correlation effect when the g-vector 
is described as g{k) - C{Vx{k),Vy{k),Q'\. This antisymmetric 
spin-orbit coupling is allowed in crystals having D„, Ci, or 
C 2 point group symmetry, such as the non-centrosymmetric 
superconductor Ulr."^^ Generally speaking, the SFM is not 
renormalized by the electron correlation when the g-vector 
of antisymmetric spin-orbit coupling is linearly related to the 
velocity of quasiparticles. This condition is satisfied in many 
theoretical models adopted for non-centrosymmetric systems. 
For instance, two-dimensional electron gases formed on semi¬ 
conductor heterostructuresand oxide interfaces'^’'^ have 
been studied on the basis of the model with the isotropic 
dispersion relation e{k) - k^jlm and RSOC with g(k) - 
(-ky, kx, Then, we indeed see the relation g{k) - 

m[-Vy{k), X’xik), 0]. Cold atom gases with a tunable synthetic 
spin-orbit coupling^" also satisfy the condition. One-, two-, 
and three-dimensional Rashba-Hubbard models satisfying the 
condition have been studied.On the other hand, 
crystals having a Td point group symmetry do not satisfy the 
condition because the g-vector of Dresselhaus-type spin-orbit 
coupling^' is represented by cubic terms with respect to the 
momentum near the F point of the Brillouin zone. 

The relation of importance, g{k) - C[-Vy(k),Vxik),Q], is 
also broken in orbitally degenerate systems. According to the 
derivation of RSOC based on multiorbital models,^^'^^’"*^’^^ 
the relation is approximately satisfied when the orbital degen¬ 
eracy is substantially lifted by a large crystal electric field. 
Otherwise, both the quasiparticle velocity and the g-vector 
of RSOC acquire a complicated momentum dependence, and 
thus they do not show a linear relation.For example, a 
cubic RSOC has been observed in the two-dimensional elec¬ 
tron gas on a SrTiOa surface,^^ and it may be induced by the 
orbital degree of freedom in t 2 g electrons.^^ Then, the ampli¬ 
tude of the SFM as well as the spin texture in the momentum 
space should be affected by the electron correlation. Thus, 
analyses of the multiorbital Rashba-Hubbard model and pe- 
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riodic Rashba-Anderson model are of particular interest, and 
therefore, we will study them in the near future. It is expected 
that the spin-orbit coupling will be renormalized through the 
renormalization of the crystal electric field, as seen in a GW 
calculation based on density functional theory. 

Finally, we comment on locally non-centrosymmetric met¬ 
als, which have global inversion symmetry but lack a lo¬ 
cal inversion symmetry on atoms. Instead of the uniform 
antisymmetric spin-orbit coupling discussed in this paper, a 
staggered antisymmetric spin-orbit coupling gives rise to a 
magneto-electric effect^* and exotic superconductivity.^^’^® 
It is expected that the staggered spin-orbit coupling will be 
enhanced by the electron correlation, as we found for non- 
centrosymmetric metals. Since the contribution of the stag¬ 
gered spin-orbit coupling is determined by comparison with 
the inter-sublattice hopping,®^ it is an important future issue 
to calculate their renormalization. 

After we submitted the first manuscript, we became aware 
of a theoretical work®^ on a related subject. Reference 58 
showed a substantial deformation of spin-split Fermi surfaces 
and the SFM near the magnetic quantum critical point, in 
sharp contrast to our results. We believe that the discrepancy 
is (at least partly) owing to the fact that the k-mass renor¬ 
malization, which plays an essential role in our conclusion, is 
neglected in Ref. 58. 
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Appendix: Fermi velocity and RSOC renormalized by 
electron correlation 


We here derive the Eermi velocity and RSOC renormalized 
by the electron correlation effects in second-order perturba¬ 
tion theory. We assume Q'|g(k)| <sc sp, as realized in most 
non-centrosymmetric metals. Then, the lowest-order theory 
with respect to the coupling constant of RSOC is justified. 

The Eermi velocity renormalized by the k-mass v^{k) is ob¬ 
tained in the zeroth order of a as 


v\k) 


ds{k) 5ReE;^(k) 


dk 


dk 


■ Aa 


d \g'(k)\ 

dk 


de(k) gRe£g(k) 
dk dk 


(AT) 


a=0 


At O' = 0, the diagonal self-energy is obtained by the Eeyn- 
man diagram in Eig. 2(2A), and thus Eo(k) = l}^^\k). There¬ 
fore, we obtain the Eermi velocity renormalized by the k-mass 
as. 


vHk) 


ds{k) 

dk 




a=0 


(A-2) 


The X- and y-components of Eq. (A-2) are expressed by 
Eqs. (18) and (19), respectively. 

The renormalization of RSOC is obtained by calculating 
the renormalized g-vector g'{k), which we obtain in the limit 


a ■ 


0 as 


g'(k) = g(k) + 


Ret^ik) 


a 


g(k)- 


dRel.'^ik) 


da 


^=0 


. (A-3) 


The spin-dependent self-energy is obtained up to the first 
order of a. 


l.{k) = 


(A-4) 


Thus, the renormalized g-vector is obtained by calculating the 
off-diagonal self-energy represented by the Eeynman diagram 
in Eig. 2(2B). Differentiating the self-energy with respect to 
a, we obtain Eqs. (21) and (22) from Eq. (A-3). 
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